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Casimir friction is analyzed for a pair of dielectric particles in relative motion. Wc first 
adopt a microscopic model for harmonically oscillating particles at finite temperature 
T moving non-relativistically with constant velocity. We use a statistical-mechanical de- 
scription where time-dependent correlations are involved. This description is physical 
and direct, and, in spite of its simplicity, is able to elucidate the essentials of the prob- 
lem. This treatment elaborates upon, and extends, an earlier theory of ours back in 1992. 
The energy change AE turns out to be finite in general, corresponding to a finite friction 
force. In the limit of zero temperature the formalism yields, however, AE — > 0, this 
being due to our assumption about constant velocity, meaning slowly varying coupling. 
For couplings varying more rapidly, there will also be a finite friction force at T = 0. 
As second part of our work, we consider the friction problem using time-dependent per- 
turbation theory. The dissipation, basically a second order effect, is obtainable with the 
use of first order theory, the reason being the absence of cross terms due to uncorrelated 
phases of eigenstatcs. The third part of the present paper is to demonstrate explicitly 
the equivalence of our results with those recently obtained by Barton (2010); this being 
not a trivial task since the formal results are seemingly quite different from each other. 

Keywords: Casimir friction; Casimir effect. 

PACS numbers: 05.40.-a, 05.20.-y, 34.20.Gj 



1. Introduction and background 

Consider the arrangement sketched in Fig. 1, namely two parallel semi- infinite di- 
electric nonmagnetic plates at micron or semi-micron separation, the upper plate 
kept at rest, the lower moving with a constant nonrelativistic velocity V . For sim- 
plicity we assume the media to be of equal composition, their permittivity being 
real, equal to e. In principle, we may allow for dispersion. The electromagnetic force 
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density in the boundary layer of the resting plate is 

f = - l -evE 2 Ve. (1) 

The essential property of the force ([lj in our present context is that it is vertically 
directed, i.e. transverse to the relative motion of the plates. And it is precisely 
this force which is able to explain all the conventional experiments in classical 
electrodynamics where dielectric media, especially liquids, are involved. Let us for 
illustration mention three examples: 

• The first is the classic radiation pressure experiment of Ashkin and Dziedzic 
D. Focused light was sent from above towards an air- water surface, and an outward 
bulge of the surface was observed, of the order of 1 fim. The light source was 
a pulsed frequency doubled YAG laser pulse (duration of each pulse 60 ns, peak 
power 3 kW, beam waist 4.5 /im). The reason for the smallness elevation is the 
large surface tension for air- water. Theoretical papers on this experiment can be 
found in Refs. [2] and [3] 

• Our second example is the experimental work of Zhang and Chang ™, showing 
clearly how a water droplet becomes deformed when hit by a laser pulse. Typical 
laser energies were 100 mJ in each pulse. Theoretical treatments for this case are 
given in Refs. [5] and [6] 

• A newer series of experiments is provided by the Bordeaux group EISI2I working 
with a two-fluid system of surfactant-coated nanodroplets in oil microemulsions near 
the critical point, implying that the surface tension can be made about 10 6 times 
smaller than the usual air-water tension. The surface displacement can accordingly 
be very large, about 70 fxm.. 

In our context, the point to be emphasized is that in all these cases the ex- 
pression ([1]) provides the correct electromagnetic force. And this brings us to the 
present theme: how can the Casimir friction force be accounted for, by the use of 
conventional electrodynamics? Obviously, the transverse force |lj cannot be of im- 
portance for friction. It is clear that we here have to leave the assumption of a real 
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permittivity. Dissipation has to be rooted in energy dissipation, that means, in a 
permittivity that is complex valued. 

It is worthwhile to show how one can come to the same conclusion via an alter- 
native approach, namely with the use of field theory. Let us start from Schwinger's 
source theory, giving how the electric field compone nts E j are related to the polar- 
ization source P k via a generalized susceptibility F^ ^ 

Ek(x) = J d 4 xT lk (x,x')P k (x'). (2) 

Stationarity means that we can take the difference r = t — t' as time variable. 
Causality means that t' < t. We take the Fourier transform to obtain the spectral 
representation Tik(r, r', uj). The Kubo formula relates V to the commutator: 

L lfc (r,r» =i f dTe iuT ([Ei(x),E k ( X ')}). (3) 
Jo 

The generalized susceptibility is the same as the retarded Green function: 
Tik(x, x') — Gf k {x,x'). Now take the Fourier transform of the two-point function 

(Ei{x)E k {x!)): 

(E i (r,Lj)E k (r',Lj')) = 2Tc{E i (r)E k (r')) u S(u) + u'), (4) 

where (Ei(r)Ek(r')) u , the spectral correlation tensor, relates to the retarded Green 
function as 

( J B i (r) J B fe (r / )) w =ImG^(r,r',u;)cothQ/3^, (3 = l/(k B T) (5) 

(this is the fluctuation-dissipation theorem). Unless the Green function has an imagi- 
nary part, the two-point function is simply zero. We thus see, also via this alternative 
route, that absorption is necessary. 

Different strategies have been followed in order to deal with the Casimir friction 
problem. The one followed here will be to adopt a microscopic model, in which we 
start from the statistical mechanics for harmonically oscillating particles at finite 
temperature moving nonrelativistically with constant velocity. Such a statistically- 
mechanical model is of course much more simple than the dielectric model shown 
in Fig. 1. We argue, however, that such a microscopic model is able to elucidate the 
essentials of the Casimir friction. 

The presentation presented here is based on our three re cent p apers I H | 12 | 13 l j n 
turn, they are based on two papers of ours some years ago EE] microscopic 
model ha s been analyzed by other investigators also, in particular by Barton recently 
| 16 | 17 | 18 l ghajj have the opportunity to compare with some of his results below. 

The other main avenue of approach in order study the friction problem, is to 
adopt the macroscopic model with permittivity properties. Let us finally refer to 
some works in that direction: Refs. fT9 l l20 l [2T | [22 } [23 l [24 } l25l 
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2. Friction force: Statistical mechanical method 

Consider a quantum mechanical two-oscillator system whose reference state is the 
one of uncoupled motion corresponding to the Hamiltonian Hq . The equilibrium sit- 
uation becomes perturbed by a time dependent term —Aq(t), where A is a time in- 
dependent operator and q(t) a classical function of time whose explicit form depends 
on the specific properties of the system. The Hamiltonian becomes H = Hq — Aq(t). 
We put 

- Aq{t) = i){r{t))x x x 2 , (6) 

where ip(r) is the coupling strength (i.e. ip is the classical potential between the 
oscillators). The separation between the oscillators is r, and Xi,X2 are the internal 
vibrational coordinates of the oscillators. When the oscillators move with respect to 
each other the coupling has to vary in time. With nonrelativistic constant relative 
velocity v the interaction will vary as 

- Aq(t) = [ip(r ) + V^(r ) • vt + ...]x lX2 , (7) 

when expanded around the initial position r = ro at t = 0. The force between the 
oscillators, called B, is 

B = -(Vip(r))xix 2 . (8) 

Note that in a mathematical sense the expansion ((TJ) requires vt to be small. Phys- 
ically, we assume nevertheless Eq. ([7]) to hold for all times, so that the interaction 
energy is proportional to t for all values of t. The natural opportunity of choosing 
q(t) — t in the interaction ([6]) thus has to be modified: as will be shown below, a 
convergence factor will be needed. 

Another point worth noticing is that the expression ([6]) corresponds to first quan- 
tization only. Quantum electrodynamic processes such as emission and absorption 
of photons (second quantization) are not accounted for by the present model. They 
were considered, however, in Ref. 1151 

The equilibrium situation with both oscillators at rest is represented by the first 
term in ([7])- It gives rise to a (reversible) equilibrium force. Thus the friction must 
be connected with the second term. To simplify, we will for the moment neglect 
the first term. By this the two oscillators will be fully uncorrelated in their relative 
position r = ro. The friction force, due to the time dependence of the interaction 
([?])■ will be a small perturbation upon the equilibrium situation. This interaction 
leads to a response A(B(i)) in the thermal average of B. And this is where the 
Kubo formula comes in: 

A(B(t)) - / cj> BA (t-t')q(t')dt' (9) 

J —OO 

(note that here 4>ba means a vector), where the response function is given by 

<t>BA{t) = ^Tr{p[A,B(t)]}. (10) 
in 
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Here p is the density matrix and B(t) is the Heisenberg operator B(i) = 

e itH/h Be -ttH/h^ where B nke A 

is time independent. Now with ([7]) and ([8]), and 
with q(t) = t, expression (| 10[) can be rewritten as (f>BA(t) = G(f>(t), with 

G = (W)(v-W), 0(t) = T>{pC(t)}, 

C(t) = ^r[xxX 2 ,X X {t)x 2 {t)]. (11) 

Thus with Eq. ([9|) and g(i') = t' the force can be written as 

F = A(B(t)) = G / 0(t - = F r + F/, 

i/ CO 

where 



OO 



F r = Gt / (12) 

is part of the reversible force by which the part of the force that represents friction 
is 

.• * 

F f =-G <j>(u)udu. (13) 
Jo 

Here the new variable u — t — t' has been introduced. The F r can be interpreted as a 
reversible force since it depends only upon position. This interpretation is consistent 
with the result obtained for the dissipation below; the F r will not contribute to the 
net total dissipation. 

If one again includes tp(ro) one has 

Gt = (W#(r +vi), (14) 

where r = ro + vi is the position at time t. By contrast, expression p^|) changes sign 
when the velocity v changes sign, and it thus represents a friction force. (Observe 
that the velocity in (| 14[) merely represents the shift in position.) Equation (fT3|) is 
the same as result (2.11) in Ref. I14l4 and the Fourier transformed version of it and 



(15) 



Eq. (jl3|) above as well, is 



iG- 



duj 



where 4>{u) = J °° 4>{t)e- luJt dt (with <j>(t) = for t < 0). 

In Ref. 14 the Fourier transformed version (|T5|) was used to obtain the explicit 
expression for the friction force. Here, we instead will use a different approach based 
on the expression (|13[) . As in the reference mentioned we then need the commutator 
(|lip . This entity again follows from the properties of quantized harmonic oscillators. 

a There is a missing minus sign in Eq. (2.11) and the F r was not taken into account. 
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We introduce annihilation and creation operators a and with commutation rela- 
tions [a,,a|] = 1 (i = 1,2; other commutators vanish). As usual, a,j(t) = aje~ %u}]t 
and ay (t) = a^e luJit . With this the coordinates are 

1/2 



■r. 



2nijUj, 



(ai + oj) (16) 



where m, and (i = 1,2) are the mass and eigenfrequency of each oscillator. To 
obtain <f>(t) from (jTTJ) we first have to calculate <p(t) = ((niri 2 |C(i)|riir£ 2 )), where 
\ri\n2) = |^2) represents eigenstates with oscillators excited to levels n\ and n 2 . 
Some calculation yields 

4>{t) = D [(2<m) + 1) cos(wit) sin(w 2 £) + (2<n 2 ) + 1) cos(w 2 i) sin(wii)] , (17) 

where D = S/(2mim2WiW2)- The energy levels = (n + ^)huj. 

In Ref. [14] the expression (17]) was Fourier transformed to obtain the friction 
force as given by (|15|) . As an alter nativ e we will here use expression (|17|) directly in 
Eq. (13]). Then we get the integral ^ 

te~ vt cos^ii) sin(u; 2 i)<ii 







^ — ^(0,), 1-0. (18) 



(772 + ft 2 ) 2 (^ 2 +fi2)2 2 ft 2 

Here ill = uji + w 2 and fi 2 = lj\ — w 2 . As mentioned above a convergence factor 
e _?? * is needed, and the limit 77 — >• is taken. Then the J7 2 — term becomes a delta 
function with prefactor determined by the integral f]x 2 (rj 2 + x 2 )~ 2 dx = tt/2. 
From (17]) we also get this integral with lo\ and w 2 interchanged. This will then give 
the result (|18|) with opposite sign with respect to the f2 2 -term. Adding up we get 
the difference of the prefactors 

,1 ,1 sinh(i,3M} 2 ) 

coth(-^i)-coth(-/3fiw 2 ' 



T ' V 2 K ; sinh(i/3fiwi)sinh(i/3fiw 2 ) 

, 77 — > 0. 

3inh(i/3fiwi) sinh(i/3fiw 2 ) ' 



(19) 



Multiplying (TSJ) with (H)]) and including the factors D and G the friction force 
becomes 

# 2 (7^)(vV^) 
F / = -- 2 ■ i,2/i a* x ^i-^a), (20) 

which is also the result (3.14) of Ref. Q3] Again one notes that there is friction only 
when the oscillators have the same frequency, and j3 should be finite, i.e. T > 0. 

When (3 — > 00, the expression (2"U]) vanishes. According to the present oscillator 
model there is thus no friction force at zero temperature. An objection against this 
result may be that it is somewhat singular due to the presence of the <5-function. 
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Thus its physical significance may not be obvious. However, rj can be kept finite. This 
will smooth out the (5-function, and the fii-term in Eq. (|18[) will give a contribution 
too. Note that this will not change our conclusions about a finite friction force 
for T > 0. (For finite 77, i.e. interaction like a short pulse, there will also be a 
contribution for T — due to the fii-term in Eq. (fi"8l) .) But here we will assume r\ 
small by which the T = contribution can be disregarded. 

In Ref. [14] the result ([20]) for the friction force was derived also by two other 
methods. These methods utilized the path integral formalism of quantum systems 
at thermal equilibrium 1 ^ 11 . The path integral can be identified with a classical poly- 
mer problem where imaginary time is a fourth dimension of length /3. Thus the 
polymers stretch out in the fourth dimension and form closed loops of periodicity 
p. For harmonic oscillators the correlation function along the polymers is obtained 
in a straightforward way. With one of the methods the convolution of the correla- 
tion functions of both oscillators were needed. The resulting Fourier transform of 
this convolution was then identified with the response function <j)(ui) used in the 
expression (fT5f ^3. 

With the other method full thermal equilibrium was utilized. Then the relative 
motion of the oscillators was regarded as a harmonic oscillator motion with low 
frequency ujq — > 0. Again with the path integral one can obtain the Fourier transform 
of the response function for the relative motion. The damping of the relative motion, 
that can be related to this response function, gives the friction force, and again 
the result (f2T)]) was recovered. Thus the three methods used in Ref. Q31 as well as 
the modification considered in the present paper, all lead to the same result, in 
contradiction to some other results in the literature, for instance that of Ref. [231 

The result for the friction force was als o ex tended to the situation with time- 
dependent or non-instantaneous interaction Then the full thermal equilibrium 
method was applicable to generalize the result. With the latter interaction there 
was also a friction from the self-interaction of a moving oscillator with itself. 

Consider now the dissipation of energy. It is identified with the work done, 
during a finite time interval starting from t = with maximum velocity v when the 
position is r = ro. As t — > 00 the motion is required to die out. To accomplish this 
we again introduce the convergence factor e~ nt (77 — > 0), whereby t is to be replaced 
with q(t) = te^^. The velocity decays exponentially, 

v^vi(t) = vq(t) = v(l-r)t)e- r ' t . (21) 

For rjt > , Vi(i) will now replace v in expression (|20[) for the friction force. 
Altogether, the total energy dissipated will be 

/OO /"OO -1 

vi(t) ■ F f q(t)dt = -v • F / / [q{t)fdt = - — v ■ F f , (22) 
-00 Jo 47 7 

where Ff is given by Eq. 00 Note that the reversible part of the force F r cx t — >• 

b The minus signs are missing in Eq. (21) of Ref. 1111 
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q(t) as given by Eq. ([T3")) will not contribute to the dissipation since / °° q(t)q{t)dt 
0. 



3. Energy dissipation calculated from first order perturbation 
theory 

Our intention now is to calculate the change AE in energy by means of quantum 
mechanical perturbation theory; this was considered also in Ref. 1121 It turns out 
that the change in energy occurs to second order in the perturbation. Nevertheless, 
time-dependent perturbation theory to the first order is sufficient to find this second 
order effect. This is because the phases of the perturbed change in amplitudes, and 
the initial amplitudes of the eigenstates, are uncorrelated at thermal equilibrium. 
Thus change in amplitudes of eigenstates will be the square of perturbed amplitudes; 
i.e. there are no cross-terms. We find that AE is positive, corresponding to a friction 
force. Doubts occasionally raised in the literature about the very existence of the 
Casimir friction effect 1^51 are thus from this standpoint laid at rest. 

Making use of the expression for AE we compare the formalism of the present 
section with that of Ref. 11 and Sec. 2 above. There, the linear response via the 
Kubo formalism was used H4U5I ^ calculate the force which in turn could be divided 



into a reversible and an irreversible part. It is the latter part that is associated 
with dissipation. A satisfactory feature is that the derivation in the present section, 
although being quite different from that of Ref. [TTJ leads to the same physical result. 

To fix the notation, we start with perturbation theory for a system at thermal 
equilibrium. The wave function can be written as tp = a n tj) n , where tp n = ipn{x) 
are the eigenstates. For simplicity we here let x represent all the coordinates of 
the system. If ip is normalized, J ip*ipdx = 1, then \a n \ 2 is the probability for the 
system to be in eigenstate n. At thermal equilibrium this probability is given by the 
Boltzmann factor P n — \a n \ 2 = Z^ 1 e~P n , where E n is the energy eigenvalue of the 
state and Z is the partition function Z = e~ /3En . Let now the Hamiltonian be 
perturbed by the time-dependent interaction V(t) = —Aq(t), where A is a quantum 
mechanical operator while q(t) is a scalar function. The A is time independent. 

Due to the perturbation the coefficients a n will change. If the system starts in a 
state m there are transitions to other states given by a change in a n , Aa n = b nm . 
The b nm is given by the standard expression 

b nm = 4 / V nm (r)e iu ™ T dr, (23) 

where 

Vnmir) = J V(T)ip m dx = -A nm q(r), 

A nm = (n\A\m). Here w„ m = ui„ - w m , with u„ = E n /h. 

As mentioned above, we will assume that the perturbation vanishes after some 
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time. Then we will obtain the total change in Aa n with 

1 f°° 
b nm = --rA nm q{-uj mn ), q(u))= / q(t)e lult dt, (24) 

™ J-oo 

where the hat denotes Fourier transform. 

From a general perspective, the system may start in a combination of eigenstates 
with transitions from several states. With this, Aa n — > ^2 m ^ n a m b nm . Now, the 
state n does not only receive contributions, but gives away contributions to other 
states also. The latter must follow from the corresponding increase of probabilities 
for the other states. Omitting the latter for the moment, the perturbed coefficients 
are 

a% n = a n + Aa n = a n + ^ a m b nm . (25) 

The a n will have complex phase factors, and in thermal equilibrium one must assume 
the phases of a n and a m (m ^ n) to be uncorrected. Thus by thermal average, 
(a* a m ) = 0. With this the new probability of the state n becomes 

Pin = ( a ln a ln) = |a«.| 2 + J]] \a<m\ 2 B nm , 

m^n 

where B nm = b nm b* nm = \b nm \ 2 . The last term is the increase in probability from the 
other states. Likewise, the state n must obey a similar loss of probability to other 
states to conserve probability. The loss to other states is thus J2 m ^n \ a n\ 2 B mn . 
With Eq. (|23|) we have b mn = 6* m , by which B mn = B nm . The latter equation 
reflects that the transition probabilities between each pair of states are the same in 
either direction. With this, the resulting perturbed probability of state n becomes 

Pm = Pn + J2 ( P ™ - P n)B nm • (26) 

m 

The change in energy can now be evaluated as 

AE = ^ E n {P m — P n )B nm = ^(£n — E m )P m B nm . (27) 

nm nrn 

Utilizing the symmetry with respect to n and m we find 

A J B=|^ e-i^ En+E ^A nm sinh(i/3A m „)i?„ m , (28) 

nm 

with A nm = E n — E m , and where 

B nm = {^/^ 2 )A nm A* nm q{-uj nm )q{uj nm ). (29) 

Note that AE > 0. The dissipation occurs to second order in the perturbation. 
To first order there is no dissipation; the changes are adiabatic. 
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4. Energy dissipation from friction force 

According to the statistical mechanical approach above - Sec. 2 and Ref. 1111- the 
energy dissipation can be written in the form 

/oo poo r />t 

v(t)F f dt = - / q{t)4>AA{t - t')q(t')dt' dt; (30) 

-oo J —oo \_J — oo 

cf. Eq. (27) in Ref. [IT] (the minus sign in front of the integrals is missing in that 
reference). The quantity q(t) is most naturally connected with the position, x{t) = 
q(t), but it can also be interpreted as a a position in a more generalized sense as 
discussed in Refs. [TT1 and [P2l so that the result (f30|) has a broader applicability. We 
will now show that this is actually the case, by showing that the result (|30p is the 
same as (|28[) , obtained by means of time-dependent perturbation theory. 
With wave function representation we first have 

e-^ B ^^ n {x)e-^*( Xl ), (31) 

n 

P AA(t) = ± J2 J M^e^rn^A^ix^e^rrn^A 

nmk 

x Mx2)e- iukt ^* k (x3)dx 1 dx2. (32) 

Thus we obtain 

Tr(pAA(t)) = \Y. A nm e^ 1 A mn e~^\ (33) 

nm 

as J ipl(x)ip n (x)dx = 5kn (x3 — x\ — x), and A nm — (n\A\m). Likewise we calculate 
Tr(pA(t)A) by exchange of uj n and uj m in Eq. (|33|) . The response function becomes 

cj> AA (t) = ±Tr {p[A,A(t)]} = i^M nm ( e -^-* - e*"™*), (34) 



ih ih 

nm 



ith 



M nm - -L e -kP^n+E m ) sinh ^A nm )A nm A* nm (35) 

(recall that A nm = E n - E m = hw nm , A mn = A* nm ). The expression for M nm 
follows if one first exchanges n and m in Eq. (j33]), then adds the resulting term to 
it and divides by 2. Some manipulation then yields 

AE = i M nm w$(u)q(-u). (36) 

nm 

With ui = uj nm = A nm /h and M nm given by the expression (|35|) this is nothing 
but the result (|28p , but now obtained by time-dependent perturbation theory. Thus 
we have been able to derive the same expression for the dissipated energy in two 
independent ways. 



November 22, 2011 1:30 WSPC/INSTRUCTION FILE 

Brevik-casimirfriction-benasque2011 



Instructions for Typing Manuscripts (Paper's Title) 11 



The result (|28p (or (1361) ) may be applied to the pair of interacting harmonic 
oscillators considered in Sec. 2, but we abstain from further details here; the reader 
is referred to Ref. [12j We give, however, the final expression for AE: 



ZA_C/ = 

877 ; 

where 



877 sum (^pwi) 



1=[-Dh\ (vVV), (38) 

with D given below Eq. (IT71) . This is the same as the result (|2"2"|) . with expression 
(1201 inserted. 



We may summarize our developments of this section: 

1. The total energy dissipation was calculated for a system perturbed by a 
varying interaction. The change in energy is basically a second order effect but was 
calculated with the use of standard time-dependent perturbation theory to first 
order only, the reason being the absence of cross terms due to uncorrelated phases 
of eigenstates. The energy change was found to be positive or zero. The result agrees 
with our previous results of Ref. Ill) obtained in a different and independent way. 

2. We have in general assumed finite temperature, and initial thermal equi- 
librium. Moreover, we have assumed low velocities and nonrelativistic mechanics, 
where by photons are not present. Photons were introduced, however, in our earlier 
study [Si. 

3. The energy change AE is finite in general. This corresponds to a finite friction 
force. In the limit T — > our formalism gives, however, that AE — > for the model 
considered in Sec. 2. This result is due to our assumption about constant velocity, 
involving slowly varying coupling. For couplings varying more rapidly, there will 
also be a friction force at T = 0, due to transitions to excited states. 



5. Equivalence between different formulations 

As mentioned, there are different formulations in the literature on how to deal 
with Casimir friction. In the present section, assuming T = 0, we wish to compare 
our formalism above with that recently given by Barton Mill ( c f. a i so R e f. fjg] 
dealing with finite temperature). In particular, Ref. [TBI dealt with a two-oscillator 
model, thus in essence the same microscopic model as ours. Barton analyzed the 
system using quantum mechanical perturbation theory. The striking point is that 
the expressions he obtained are seemingly quite different from those we obtained 
in Ref. [TT] In Barton's own words (Ref. [TBJ Sect. 3) "..in view of the manifold 
current controversies about quantum-governed frictional force generally, it seems 
well worth exploring whether such differences reflect substantiate disagreement or 
only a confusion of terms". Our present investigation is a follow-up of Barton's 
suggestion. We intend to demonstrate explicitly that the obtained expressions for 
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the dissipated Casimir energy are in fact in agreement with each other, thus a 
reassuring result. 

Assume that the oscillators have the same eigenfrequency u> and the same mass 
m. They interact via a time-dependent Hamiltonian 

e 2 

Hint = ^yiy2, (39) 

(Gaussian units assumed). Here e is the elementary charge, y\ and j/2 the oscil- 
lator coordinates, and s = s(i) is the vectorial distance between the mass centers, 
varying with time because of the relative motion of the oscillators. Introducing new 
coordinates 

V± = ~ 7r , (40) 
one can write the interaction Hamiltonian as 

1 e 2 
H iDt ±=±-qy±, 9=^3- (41) 

By use of time-dependent perturbation theory the total energy dissipated is then 
found to be 

AE = 2x 2haj\c(oo)\ 2 , 

c{t) = ~k / oo ^( 2 +ly+l°+> eW ' (42) 

as given by Eqs. (3.3) and (2.4) respectively, in Ref. H6l At T = 0, only excitations 
from the ground state are possible. 

To compare Barton's result (|42j) with ours, we first have to evaluate the matrix 
elements in (|42j) . In terms of the creation and annihilation operators we have (u>± — > 
u! for small perturbations) 

y±=Vb(a ± + al), b=^. (43) 

Then, 

(2 + |y 2 |0+) = b (2 + | a ; 2 |0+) = V2b. (44) 
Together with Eq. (|4"2")l this gives 

AE = 8hujb 2 \I (oo)\ 2 , 



i 



t 

I „ 2iu>t' 



/ (*) = -^/ dt'qe^. (45) 
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Proceed now to compare this result with those that we derived in Refs. fill and [T2l 
at T = 0. The interaction Hamiltonian, written as H- lnt = —Aq(t), now corresponds 
to 

e 2 

A = -ym, and q{t)=q=—. (46) 

s 6 

Recall that A is a time- independent operator accounting for the quantum mechan- 
ical properties of the system, while q(t) is a classical function of time. At T = the 
system is in its ground state with probability equal to one. 

The change in energy is given by Eq. (|27p. where E n is the energy in the (un- 
perturbed) eigenstate n, and P n — \a n \ 2 is the probability of the system to be in 
this state. As mentioned, we start from the ground state so that P m —> Poo = 1- 
Further, as B nm = |6 nm | with b nm the transition coefficient between states m and 
n, Eq. ([27)) reduces to 

AE = (E u - E m )B 1W0 , (47) 

where En — Eqo = 2haj is the energy difference between the state |11) where both 
oscillators are excited to the first level, and the ground state |00). The coefficient 
-Bnoo is the transition probability between these two states. 
What remains is to calculate -Bnoo- Using Eq. (|24|) we get 

q(-Unm) g(-wnoo) = g(-2w) = 2ihI(oo), (48) 
with I(oo) given by Eq. (JTTJ). Further, 

A nm A U oo = (H| - yiy2|00) 

= -<1|2/i|0><1|2/ 2 |0) = -(l| yi |0) 2 = -6<1|«I|0) 2 = -b. (49) 
Altogether, when inserted into Eq. (j47|) we obtain 

Bnoo = -^\A 1W0 \ 2 q{-2u)q(2Lj) = 46 2 |/( TO )| 2 , 

AE = Shu;b 2 \I (oo)\ 2 , (50) 



which coincides with the result obtained by Barton, Eq. (|42[) above. 

As an additional remark we note that the situation with zero friction at T — for 
slowly varying forces can be analyzed in a straightforward way from the equations 
above. With time-dependent part (q = 0, t < 0), 

q(t) = te^, t > 0, (51) 

we namely find when r\ — > 



1 



?M = TZ , ( 52 ) 



1 7T 

q(uj)q(— uj) = -—^ 7T-r — > - — kS(uj). (53) 

v J y ' (f] 2 +cu 2 ) 2 2t]uj 2 y ' v ' 
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Thus with w/0, Eq. ([5B"f will give zero for the dissipated energy AE. 

[By looking at this in some more detail it would seem that the middle term in 
Eq. (|53[) . before taking the limit rj — > 0, implies there to be some dissipation. This is 
physically an artefact, due to our assumption of an abrupt change of q(t) at t = 0.] 

The equivalence between Barton's results and ours is therewith shown. We thus 
hope to have shed light on one of the subtle issues in the Casimir friction world 
demonstrating that quite different approaches can lead to the same result. 
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